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ABSTRACT 

In the Fornax dwarf spheroidal galaxy the globular clusters are distributed 
widely, without any significant central concentration. Oh et al. (2000) pointed 
out that such a distribution is paradoxical: dynamical friction effects estimated 
using single-component King models would have forced the globular clusters to 
spiral down to the center of the galaxy well within a Hubble time. This paper is 
devoted to a discussion of this paradox. We describe a model in which the stars 
of the dwarf spheroidal galaxy are embedded in a cloud of dark matter, and each 
of these components is specified by its own phase space distribution function. 
This model allows us to fit self-consistently the observed luminosity profile and 
the spatial variation of the velocity dispersion of the stars. This fitting yields 
two basic parameters, related to the central density and velocity dispersion, that 
characterize the phase space distribution of dark matter. The dynamical friction 
effects calculated on the basis of this self-consistent model are small enough that 
the observed spatial distribution of the globular clusters poses no difficulty, and 
the apparent paradox is resolved. Thus we have at hand a model for Fornax that 
reproduces the main observed features of this dwarf spheroidal galaxy. 

Subject headings: dark matter - galaxies: dwarf - galaxies: kinematics and dy- 
namics 

1. Introduction 

The masses and the phase space structure of the smallest stellar systems that form 
in the Universe provide the most direct clues to decipher the nature of dark matter and 
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the process of galaxy formation. Among these systems, dwarf spheroidal galaxies are the 
most suitable for the study of dark matter. They have low densities of visible matter in the 
form of stars (about a million times smaller than those encountered in globular clusters of 
comparable masses), and both their internal and external dynamics are dominated by the 
cloud of dark matter (more commonly called the halo) in which they are embedded. Their 
stabihty against tidal disruption and also their internal dynamics depend upon details of the 
structure and extent of the dark-matter cloud, which in turn is sensitive to the phase-space 
structure of the dark matter particles constituting the halo. 

In recent years, improved astronomical observations of these faint systems have become 
avail able for comparison with theoretical studies. One of the earliest studies of these systems 
is by lFaber &: LinI (119831 ) . who used luminosity profile observations to show that dwarf galax- 
ies are dominated by dark matter even in their cores, drawing attention to dwarf sph e roida l 
galaxies as candidates for the study of dark matter. Around the same time, lAaronsoru (119831 ) 
measured the radial velocities for several carbon stars in these systems with sufficient accu- 
racy to support the idea th at dark matter plays a domina nt role in their internal kinematics. 
Prompted by these studies, ICowsik &: GhoshI (1l986l . 119871 ) developed an "embedding model" 
where the dwarf spheroidals are embedded in an extensi ve cloud of dark matter. A similar 
approach was championed by iPrvor fc Kormendv I (11990 ) in t he context of dwarf spheroidal 
galaxies. Tremaina ( 1976 ) and Hernandez fc Gilmorel ( 1998 ) discussed dynamical friction 
effects on globular clusters in dwarf galaxies. They derived simplified analytical formulae for 
the orbital e volution of a massive body in a dark h alo which is assumed to follow the King 
distribution (jKinglll966l : iBinney fc Tremaind l2008l ). finding that the dynamical timescales 
for the globular clusters to sink to the center of the galaxy are about 1 Gyr, and that these 
timescales increase as the core radius increases. 



Of particular relevance to the present study is the work of lOh et al.l (120001 ). who pointed 
out that the five globular clusters gravitationally bound to the dwarf spheroidal galaxy in 
Fornax paradoxically "preserve their diffuse spatial distribution despite the fact that the 
clusters' orbital decay t imesc ale is much shorter than the estimated age of the host galaxy." 
Sanchez-Salcedo et al.l (120061 ) have shown that the problem is exacerbated if modified New- 
tonian dynamics is applicable. There is a sim ilar problem wit h the orbital decay tim escales 
in dwarf elliptical galaxies (for details see e.g. iLotz et al.ll200ll ). iGoerdt et al.l (120061 ) consid- 
ered the effects of dynamical friction in da rk matter halos with profiles haying a central cusp. 



as su ggested by cosmological simulations (lHernquistlll990l : iNavarro et al.lll997l : iMoore et al. 



19991 ) and showed that "in a cuspy CDM halo the globulars would sink to the center from 
their current position within a few Gyrs, presenting a puzzle." They further point out that 



should the stellar population be embedded in a dark matter ha . 



ficulty would be resolved. In making this point, iGoerdt et al.l (120061 ) assume that the dark 



o wit h a large core, this dif- 
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matter halo has a density profile given by 



p{r) = pD{m + {r/rs)\ 



(1) 



with a central density Pd(0) = O.IMqPc^^ and scale size = 2.4kpc. They also assume that 
the density sharply cuts off at ~ 50kpc. Thereafter Jeans' equations representing the balance 
between gravitational forces and the gradient in pressure provides them with the (f^(r)) for 
this assumed density profile. This allows them to estimate the time scales for the dynamical 
friction to operate and conclude that in such a halo the dynamical friction effects on the 
globular clusters are small. This immediately opens up the following questions: (i) Does the 
density distribution for dark matter given in Eq. [1], with the chosen values for the parameters 
P£)(0) and Vg, represent correctly the density profile of the dark matter halo in Fornax? (ii) 
Would the stars embedded in such a halo follow the observed spatial distribution? (iii) 
Would the observed v^:rms{r) of these stars be correctly reproduced? 

In other words, in order to explicitly resolve the problem posed by the wide distribution 
of globular clusters in Fornax, we should work within the context of a model that will allow 
us to formally derive the density profile p*(r) of the stars and the profile of their velocity 
dispersion v^:rms{r) in terms of the parameters describing the phase space distribution of 
dark matter. Such a phase space distribution, with the best choice for the parameters, may 
then be inserted into the relevant integral in Chandrasekhar's formula to estimate the time 
constant for the m igration of globu l ar clusters toward the center o f For nax. As a foUo A V up o f 
the earlier work of iTremaind (Il976l ). [Hernandez fc Gilmord (119981 ) . and lGoerdt et al.l (120061 ) . 
this paper describes an attempt to construct such a model and address the problem posed 
by the wide distribution of globular clusters in Fornax. 

Fortunately aiding the model building effort, astronomical observations of the dwarf 
sphero i dal gal axy in Fornax have improved considerably since the pioneering studies by 
Hodgel Jl961al lbl). and now excellen t data sets are available. The galaxy is located at a 
distance of ~138 kp c (IMateol Il998[). It covers a large area of the sky, with an estimated 
tidal radius of ^71' (llrwin &: Hatzidimitrioulll995l ). corresponding to ~ 2.85 kpc. A careful 



wide- field survey by I Coleman et al.l (120051 ) provided a precise radial number density profile of 
the RGB-selected stars, showing that the observed distribution increasingly deviates from a 
single-component King distribution beyond the core radius (see Fig. [1]). The surface density 
beyond the core radius falls off as a power law and s ignificant den s ities w ere noted up to 90', 
far beyond the earlier estimates of the tidal radius. Walker et al. ( 20061 ) measured the radial 
velocities of a large number of stars in the galaxy and derived the radial velocity dispersion 
(v^(i?))^/^ as a function of the projected distance R. They noted that the radial dependence is 
fiat, again deviating significantly from the decreasing behavior as a function of R that would 
be expected for a single-component King distribution (Fig. [2]). The five globular clusters 
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that pose a problem for dynamical friction timescale s are located at pr ojected distances 
R = 0.24, 0.43, 1.05, 1.43 and 1.60 kpc, respectively JWalker et allbood ). Of course their 
actual radial distances would be larger, by a factor of ~ a/3/2 on the average. 

In order to model the observations noted above, we derive in Section [2] the radial depen- 
dence of the velocity dispersion and the surface density profile of stars within the context 
of a dynamical model, by formulating the Boltzmann-Poisson equations in the embedding 
approximation. The solution to these equations is obtained through a standard procedure, 
which is briefly described. We assume that the underlying phase space distributions, both 
for the dark m atter and fo r the stars, are nearly isothermal and follow the lowered isothermal 
distributions (lKing||l966l ). In making this choice we were guided by extensive earlier work 
adopti ng lowered isothermal phase space distributions to model these systems (see for ex- 
ample llrwin fc Hatzidimitrioul (119951 ) , iBinney fc Tremaind (120081 ) and references therein) , as 
well as by the requirement of simplicity. Since there are two components in the Fornax system 
- the dark matter particles and the stars - we need at least two distributions to specify them, 
and the reduced isothermal distributions have the advantage that they form cored systems 
with finite masses and have been adopted extensively before for modeling various astronom- 
ical systems, including the dwarf spheroidal galaxies. It is perhaps relevant to point out here 
that rece ntly some concerns have been expressed regarding their applicability to galactic 
systems ( IGilmore et al.ll2007l ). These authors have noted that the two-body relaxation time 
in these systems, T2 ~ {SN/ log N)Tcross, is much longer than the age of the universe, and 
as such the system would not have relaxed to a nearly Maxwellian phase space distribution. 
This remark is true and the absence of any segregation of massive stars into the central 
regions (except in the case of the densest globular clusters) indeed points to the absence of 
such equilibrium induced by two-body interactions. However, the smooth density profiles of 
these systems and our ability to model them with reduced isoth e rmal distribution functions 
beg for an explanation. The pioneering study by iLynden-Belll (119671 ) has shown that the 
equilibria in such systems are dictated by collective effects. He showed that even though the 
'fine-grained' distribution takes almost infinite time to evolve, the 'coarse-grained' distribu- 
tion function rapidly relaxes due to collective effects into a quasi-stationary state - a process 
called 'violent relaxation'. His work triggered wide-ranging studies of the collective effects 
of long-range forces. As long as initial conditions avoid parametric oscillations and kinetic 
energy approximately balances the potentia l energy, 2T + U p a 0, the coarse grained distribu- 
tion does turn out to be quasi-Maxwellian (ILevin et al.l 120081 ). It may be noted further that 
it is the fluctuation of the gravitational potential due to collective effects that brings about 
this quasistationary state: all stars acquire the same distribution of velocities and there is 
no segregation, in contrast to equilibria driven by two body interactions. Accordingly we 
adopt reduced isothermal distribution functions for the stars and dark matter particles, solve 
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self-consistently the equations of the embedding model and therefrom compute the slowdown 
rate due to dynamical friction. In Section [3] we compare the theoretical predictions with the 
observed velocity dispersion and luminosity profile, and we show that the model indeed fits 
the observations very well. Thus the main parameters of our model (a density parameter 
pcD and a velocity dispersion parameter an) are determined with sufficient accuracy, allow- 
ing us to estimate with confidence the time needed for the globular clusters to sink to the 
center due to dynamical friction. This turns out to be longer than the Hubble time, thereby 
resolving the paradox. Section H] is devoted to a general discussion and concluding remarks. 



2. The Dynamical Model 

We begin the theoretical analysis by developing a dynamical model for the dwarf galaxy 
in Fornax. A comparison of the predictions of the model with the observations of the profiles 
of the root mean square of the line of sight velocities Vrms*{R)i and the projected number 
density, of the stars shows that the model provides an excellent fit to the observations 

and provides the best values of the two basic parameters pen and characterizing the 
phase space distribution function of the dark matter particles in Fornax. Using these we 
analytically integrate the incomplete integral over the phase space distribution appearing in 
the Chandrasekhar formula for dynamical friction. 

Let the phase space distribution both of the dark matter fa and of the stars /* be of 
the 'reduced isothermal' form (King 1966) 



for > 



= for < . (2) 

The subscript i could either be 'D' or for the particles of dark matter or the stars as 
appropriate, and the 'binding energy' e is given by 



Si 



where 0(r) = (pnif) + 4>*{f') is the total gravitational potential deriving contributions from 
dark matter and the stars. The terminal radii rxi for the two components are indirectly 
specified by the choice of the value of the total potential at these locations. Since the total 
gravitational potential due both to the dark matter and the stars occurs in both the distri- 
bution functions, it couples the two components. Notice that these distribution functions 
depend only on f ^ which implies isotropy for the velocity distribution, and the distribution 
functions, being functions of energy alone, satisfy the stationary collisionless Boltzmann 
equation. The parameter r^i represents the apogalacticon at which the density vanishes. 



- 6 - 



and this is indirectly specified by the ch oice of [(f>{rKi) — (f){0)]. The integral over the velocities 
(lKing||l966l : iBinney fc Tremaind l2008l ) provides the spatial density poi^) and p^:{r): 



Pi = Pa 




(4) 



where yi 
equations 



^i/crf are the scaled potentials satisfying the spherically symmetric Poisson 



2 dyi ^ d'^yi 
r dr dr'^ 



AttG 

-Pi{r) ■ 



(5) 



Notice that both pnir) and p*(r) are functions of the total potential 0(r) = 4>D{r) + 0*(r). 
This represents the coupling between the two systems; as the density of one of these be- 
comes higher its contribution increases and the total potential becomes deeper, making both 
components more compact. Accordingly, these equations represent two coupled non-linear 
differential equations, which are solved by Taylor expansion of y^ about the origin, say up 
to the sixth order in r, and thence integrated using standard numerical methods (an adap- 
tive Runge-Kutta scheme) to extend the solutions from small values of r to any desired 
distance. Noting that the mass to light ratio of dwarf spheroidals, including that in Fornax, 
are very large, we are able to simplify the numerical analysis by neglecting the contribu- 
tion of the stars to the overall gravitational potential (f){r). From the numerical values of 
yD and dyn/dr all quantities of interest may be evaluated: for example the circular ve- 
locity Vc{r) = {ajjr\dy£,/dr\}^^'^, and the mass of dark matter contained within a radius 
Mi:){r) = a'jjr'^\dy D / dr\ / G . Furthermore, noting that with dominance of dark matter 



y*{r) = [yoir) - yoirK*)] (^1/(7 



2 



(6) 



the density of the stars is evaluated using Eq. |H Now, the projection of any variable Z{r 
on the plane of the sky is given by 

rdr 



(7) 



where R is the projected radial distance. With this, the surface density profile of stars is 
given by 

K{R) = S{p,{r),R). (8) 
Defining the second moment of the phase space distribution of stars as 

/2er 



TJr) 



3pc*cr* 



e^*erf(yy2) _ 



5/2 



3/2 ^ 1/2 



(9) 



- 7- 



and recalling the assumption of isotropy, 

TJr 



(vlir)) 



3p*(r) 



and Vrjns,*{R) 



1 S{T,,R) 
3 S{p^:, R) 



1/2 



(10) 



In Eqns. [TlfTOl we have the theoretical functional forms for the profile of stellar surface 
density S*(_R) and the profile of stellar velocity dispersion Vrms*{R), which may be compared 
with the corresponding profiles obtained from observations to determine the best choice of 
the parameters entering the distribution functions. 



Now, Chandrasekhar's formula (j Chandr asekhan Il943l : iBinney fc Tremaind |2008| ) pro- 
vides an estimate of the dynamical friction acting on the globular clusters in Fornax: 



dV 
~dt 



V 



\V\^ 



(11) 



where we have set (M + m) ~ M, and estimate the Coulomb logarithm as In A = 4; also here 
Vm is the smaller of V and [2^(r)]^/^, the latter being the escape velocity from the system. 
Defining = K^/2a^ we get 



dV _ _V_ 



X 



A-nG'^Mp^ In A 



(12) 



e^^ <^ erf (x^) - 



TT 



3v^ 



Thus, we have at hand the formalism needed to fit the observed profiles of number 
counts and of dispersion in the stellar velocity and obtain the parameters pz)(0) and an 
crucial to the phase space distribution function of the dark matter particles. These will 
allow us to use Eq. [12] to evaluate the dynamical friction effects on the globular clusters in 
Fornax. 



3. Comparison with observations 



The theoretical profiles Vrms* {R) and S^, (i?) depend on the choice of the model parame- 
ters: pcD (or more conveniently p£)(0)), an, tkd for the dark matter and pc* (or p*(0)), 
and tk* for the stars. Amongst these p=i<(0) provides the normalization constant for S^,(/2) for 
the particular type of stars which are being counted, and tk* the location where we would 
like to cut-off the stellar distribut i on. N oting that the farthest data point in the stellar 
profile observed by IColeman et al.l (120051 ) lies at a projected distance R ~ 4kpc, we may 
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choose rx* to be any value significantly greater than this, say lOkpc. The observed profile is 
insensitive to this choice. Similarly the dependence of r^D is also weak, fixing the total mass 
of the dark matter halo. Thus the crucial parameters that control the theoretical predictions 
are a^:, an and Pi?(0). It may be appropriate to point out here that ^/^(O) and 2/*(0) are not 
additional independent parameters. The choice of Pd(0), cr^i, tkd and corresponding stellar 
parameters fixes them through the requirement that UDirKD) = and y*{rK*) = 0. 

We start by fitting the observed luminosity profile (R) based on the embedding model 
presented in the earlier sections. Our fit is shown in Fig. [T]for the following set of parameters: 

Pd{0) = O.O4M0pc"^ cxd = 19.5km ■ s'^ a* = 12km ■ s"^ . (13) 

Note that the model fits the observations very well, including the power law tail at large 
rad ii. In contrast, a sing le component King model does not fit the observations, as noted 



by I Coleman et al.l (120051 ). We emphasize that the fit to 'E^{R) is very sensitive to the 
choice of Pd(0), (Jd, and a^, and at most a ten percent variation may be accommodated 
in the values of these parameters given in Eq. [T3i The values of Po(0) and a^, determine 
approximately the core radius of the stellar distribution in the embedding model; this radius, 
ro* = [9a'^ / AttG pd{0)Y^'^ = 650pc, marks the location where the projected stellar density 
falls to about one half of its central value. It is worth noting here that the quantity (2(T^/cr^) 
yields approximately d\iap^{r)/dhi{r), i.e. the power law exponent of the fall off of p*(r) 
beyond the core; or equivalently, [(2cr^/(T^) — 1] gives the index for S^,(i?). 

In order to ensure that our model is consistent with all the rest of the available ob- 
servations we calculate (f^^(i?))^/^ with the same set of parameters fix ed by fitting S^(-R) , 



and compare the prediction with the observed radial velocity dispersion (IWalker et al.ll2006l ). 
The data and our fit are shown in Fig. [21 Even though the statistical uncertainties in the 
data set are large, our model is consistent with the observations. On the other hand, as 
pointed out by Walker et al., their observations are inconsistent with the single component 
King model (dashed line). 

Having fixed the parameters of the model, it is a straightforward matter to calculate 
dynamical friction effects on the globular clusters using Eq. [121 A convenient way to display 
the results is to assume circular orbits at various radial distances of closest approach and 
calculate the dynamical friction rate for V = Vc, where Vc is the circular velocity. There are 
two reasons for doing so: (1) if dynamical friction effects were indeed significant, then the 
orbits of the globular clusters would tend to be circularized by friction; (2) dynamical friction 
effects are maximized in the circular case, because non-circular orbits with the same distance 
of closest approach have higher velocities, and they also sample the outer, low-density regions 
of the dark matter distribution. 
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Choosing the best-fit parameters, p_d(0) = O.O4M0pc ^ and au = 19.5km ■ s ^, we 



calculate the rate of dynamical friction given in Eq. [12] and plot it in Fig. [31 Very close to 
the center (i.e., for small r) the dynamical friction formula ([T^ simplifies to 



which implies a lifetime 



T 



dV 
'dt 



V 



167r2G2^1nAp^(0) 

3(27r(T2)3/2 



V . 



(14) 



dV/dt 



4v^G'2MlnAp^(0) 



4.5 Gyr 



(15) 



for globular clusters in Fornax. This simple estimate is consistent with the exact calculations 
presented in Fig. [31 The effective lifetime, r, is ~ 4.5 Gyr in the central parts of the dwarf 
spheroidal galaxy, it is nearly constant out to the core radius, and for larger radii it increases 
roughly quadratically with R. Since it would take several e-folding times for a globular 
cluster to sink to the center, even a globular cluster initially formed in the core region would 
not have migrated to the center within a Hubble time. For globular clusters that condensed 
in the outer regions or with elongated orbits, the effects of dynamical friction are much 
weaker than this estimate. 

To be more specific, the dynamical friction causes loss of angular momentum at a rate 

dv 



dL 
~dt 



dt 



(16) 



Here the angular momentum L is taken to be r ■ Vc, and the dynamical friction dv/dt is 



evaluated using Eq. [T^l with V 



keeping in mind that friction will tend to make the 



orbits circular. For a given perigalacticon circular orbits suffer the highest rate of dynamical 
friction, as elliptical orbits will spend a large part of their orbital time in the lower density 
outer regions. In order to calculate the radial migration rate, we write 



dt \dt)'\dr) \dt)'dr^ 



(17) 



and define the radial migration time constant by 

dr 



T,{r-) ^ r/ 



dt 



(1^ 



This quantity is plotted in in Fig. [51 Note that becomes very large in the outer 
region, and becomes equal to r^, well within the core radius. Indeed, the radial migration 
rate given in Eq. [17] is easily integrated to obtain r{t) for various assumed values of r„. 
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the present location of the globular clusters in Fornax. We show in Fig. [6] the past and 
the future trajectories, starting from the time of their formation ~ 10^ million years ago, 
arriving at their current location in the present epoch, and moving into the future another 
lO'^ million years. These trajectories exclude any paradoxical behavior and are consistent 
with the astrophysical assessment of the regions of their formation and subsequent evolution 
of their orbits. They show that the globular clusters that were formed at distances between 
~ 1 and 2 kpc from the center of Fornax, suffered dynamical friction at the expected rate 
and arrived at their current locations at the present epoch. During the next 10 billion years 
they will migrate closer to the center, but will remain distinct without forming a nucleus. 

Thus the apparent paradox posed by the d iffuse distribution of globular clusters in the 
Fornax dwarf spheroidal galaxy (jOh et al.ll2000l ) is resolved. In order to provide some of the 
derived parameters of the stellar component in F ornax, we take the total visible magnitude 
of Fornax rriL ^ 9.3 (IKarachentsev et al.l |2004| ). distance d = 138kpc and a bolometric 



correction ~ 1"^ to estimate the total luminosity of L ^ 7.5 x IO^Lq. (For a bolometric 
correction of ~ 0.5™, L ^ 6 x IO^Lq). In Fig. H] we show the density of the stars p*(r) 
in arbitrary units, the luminosity within a radius r and the mass to luminosity ratio as a 
function of the radius; of course the mass is dominated by the dark matter. 



4. Discussion and conclusions 

The dynamical friction suffered by astronomical objects (such as globular clusters) de- 
pends not only on the local spatial density of the surrounding region, but also on its phase 
space structure. Since the phase-space density is invariant along a dynamical trajectory, 
this sensitivity implies that several parameters of interest (such as the dispersion velocities 
and the extent of the dark matter halo) may be determined by examining the role played 
by dynamical friction in extended structures in the Universe. These remarks are illustrated 
here through a study of the dynamics of globular clusters surrounding the dwarf spheroidal 
galaxy in Fornax. To summarize, a model based on reduced isothermal distributions both 
for the stars and the dark matter provides a very good fit to the radial velocity dispersions 
Vrms*{R) and the luminosity profile £=„(/?) observed for the stars in the dwarf spheroidal in 
Fornax. These fits determine Pd(0) and an, the two crucial parameters of the phase space 
distribution of dark matter. The rates of dynamical frictio n calculat e d usin g these parame- 



ters are sufficiently small that the paradox pointed out by lOh et al.l (j200d ) is resolved. An 
application of such models to other systems including galaxy clusters would provide further 
insight into the phase space distribution of dark matter and the process of structure forma- 
tion in the universe. It would also be useful to attempt self-consistent fits to the observations 
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of Fornax, using the phase space distributions obtained with numerical simulations, as soon 
as they become available. 
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Fig. 1 . — The theoretical ly calculated surface d ensity in Eq. |8] (solid black line) is compared 
with the observations by I Coleman et al.l (120051 ) and their single component King model fit 
(dashed line corresponding to a* ~ 11km- s~^ and r^* ^ 2.6kpc), which progressively falls 
below the observed surface densities at distances beyond the core radius. The embedding 
model (solid black line) provides a very good fit. 
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Velocity Dispersion 



R[pc] 



Fig. 2. — The observed dispersion in the radial velocities using 182 stars (IWalker et al.ll2006l ) 
is compared with our theoretical model (solid black line). A single component King model 
presented by Walker et al. (dash-dot line) with 0(O)/o"^ = 3.26 and tidal radius r^* ^ 3kpc 
predicts smaller dispersion velocities beyond the core radius. 
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Dark Matter Properties 
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Fig. 3. — The lifetime r„ = \V/{dV/dt)\ for orbital decay due to dynamical friction, in Myr 
(dash-dotted line) is calculated for circular orbits, using the best-fit parameters for the dark 
matter distribution obtained from the observed stellar velocity dispersion and stellar density 
profiles. For completeness we also show the enclosed mass in dark matter in units of 
(dotted line), the circular speed Vc in km/s (dashed line) and the density of dark matter 
in Mq/pc^ (solid line). 
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stellar Properties 
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Fig. 4. — The theoretical estimates of the stellar density p*(r) in Mq/pc^ (solid line), stellar 
luminosity L enclosed within a radius r in Lq (dotted line) and the total mass to luminosity 
ratio M/ L also enclosed within r in units of Mq/Lq (dash-dotted line) are shown as a function 
of the radial coordinate. Note that the density of stars becomes less than one thousandth 
of its central value for r > 3kpc, and it is unclear whether these stars really belong to the 
dwarf galaxy. This is indicated by the dashed-line extension to the model estimates of p*(r). 
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Radial Migration Time Constant 
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Fig. 5. — The radial migration time constant r^. = r /\dv/dt\ is displayed as a function of r. 




Fig. 6. — The radial trajectories of globular clusters from the time of the formation up to 
lOGyr in the future is shown. The trajectories pass through their present location, r = 
A/3/2i?, lOGyr after formation. 



